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DESCRIPTION OF INVENTION: 

An optimization process for roads planning. 
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ABSTRACT: 

Using optimization process, through linear or nonlinear programming, to obtain the 
optimal- points of intersection, vertical alignment of roads and lots, cross-sections 
Quantity cost estimation and work program of excavation material transportation for the 
entire roads system and subdivision, according to topography, geology roads standards 
and earthworks prices. This invention utilizes the information generated by a computer- 
aided design (CAD) program as input to the optimization algorithm. 

14 Claims, 4 Drawing Sheet 

BACKGROUND OF THE INVENTION: 

Before the present invention the design for the roads was done separately for each road in 
the system The design of a roads system usually began with the planning stage where the 
civil engineer arranges the road to satisfy the standard specifications such as grade and 
meeting of centerlines at intersections. Then a draftsman produces an initial design using 
one of numerous software programs to draft profiles, vertical alignment, and cross- 
sections calculate the cut-and-fill quantities. Usually the initial design needs to be 
rearranged numerous times by the engineer to satisfy all of the physical constraints. This 
iteration may be repeated numerous times before the road system satisfies all of the 
physical constraints and balance the cut-and-fill quantities. The design has been typically 



taken on account. The optimization approach was not be used for roads des. = n. 

OR IECTS OF THE INVENTION: 

"fobjert of this invention is to develop a process of road design that arranges optimally 
Lometr; of roads system that satisfies the roads standard applicable into the most 
economical entity. 

The object of this invention is also to develop a process to compare different layouts and 
lifferem unit prices to obtain the cost-effectiveness in the preliminary planning stage. 

nF^rRTPTION OF THE INVENTION: 

The crimen- modeling system establishes a direct relationship between the current 
the optimization algorithms (Fig. 1). 

The optimization is reached by three methods that differ each from other by number of 
variables, time of performance: 

a Grid optimization (Nonlinear programming approach) - : VP"***' s for 
^homogeneous topography, when geological data is availab ^ *^£e 
and in the roads system exist fixed points (see paragraph 5) The °P^ za *° n 
Obtained in one phase. This approach is time consuming because it handles 

many variables. 

b Centerline optimization (Linear programming approach) - aPP^P^te for 
^ homogeneousVography, when geological data is not J^^^J? 

objective function represents difference between cut and fill. The optimization 

obtained in two phases. This approach is time consuming. 

c Centerline optimization (Nonlinear programming approach) - appropriate for 

eStSg elections along centerline.. The optimization is obtamed ,n two phases. The 

result is obtained faster. 



1. Grid optimization (Nonlinear programming approach): 

The optimization is reached through nonlinear programming when objective is 
earthworks total cost. Variables are differences between proposal roads elevations and 
existing topography elevations at all grid points on right of way. Constraints are present 
roads standards. 

Grid points (GPs) are 3D points at nodes of rectangular grid of definite size (e.g. 5m x 
5m) that located on right of way. Boundary grid points (BGPs) are GPs located on edges 
of roads, i.e. points for which distance to road centerline is equal to half of road width 
(with some tolerance). BGPs are considered for retaining walls estimation. 
Control points (CPs) are 3D points on proposal centerlines of roads system (e.g. for each 
100m). CP might not be GP. (Fig. 2). Grade of road can be changed only at CPs. 
Coordinates X and Y of CPs and GPs are entered as input data and coordinates Z 
are an output of optimization procedure. Recommended locating CPs at points of 
centerlines where grades of existing profile are changed. This will improve optimization 
process. Using CPs and GPs elevations as decision variables we now define the nonlinear 
programming task: 

Let Q(X'q, Y'q,Z'q) be a GP. Therefore Q lies on some road of given roads system. Then 
there must exist A(Xa,Ya,Za) and B(Xb, Yb,Zb) - adjacent CPs on this road centerline so 
that Q is located between A and B. It means that if Q' is nearest to Q point on road 
centerline including A and B, then Q is located between A and B. Points Q and Q' must 
have same elevations. Therefore for any point Q from GPs exist X, p. e [0,1] so that: 
Z'q = X*Za + ^*Zb, 

Where 

Z'q, Za, Zb are proposal elevations of points Q, A, B respectively. 
If interval AB is straight then X = (AQ,AB)/(AB,AB) € [0,1], 

It means that each GP proposal elevation may be received via linear interpolation from 

proposal elevations of corresponding pair of CPs. 

Let 

Zq = Z'q - Ze, where Z'q is proposed elevation and Ze is existing elevation at point Q. 
Zq is relative elevation at point Q. Zq is a decision variable and used in the objective 
function. 

Objective is earthworks total cost: 

q=N 

(1.1) Z ( Fq(Zq)*G+ Wq(Zq) + Cq(Zq)*G + Vq(Zq))*G + D 
q=l 

where 

G is grid size (e.g. 5m); 
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Fq(Zq) is cos. of fill a, point Q, where proposal GP elevation is btgger than existing CP 
WqTzqTiscostrf retaining wall for fill a. point Q, where proposal GP is above existing 

uSZZl^erence between fill and excavation volumes moving in or out of 

T^<°«tTwa(Za) = Fq(Zq) = 0- also if Q is not BGP then Wq(Zq)=0; 
ff In > 0 hen CaSa) = Vq£q = 0;' also if Q is not BGP then Vq(Zq)=0. 
^d^« tottotlS comprise the overall tota! earthworks cost are now 

Lettnce of fill = fO (money unit / volume unit), then Fq(Zq) could be written so: 

Fq(Zq)=0.5*f0*(|Zql + z <0> 

Hence, if Zq < 0 then Fq(Zq) = 0. 

If excavation price is constant-cO, then 

Cq(Zq) = 0.5*c0*( | Zq | - Zq 

^Z££%2Z%*«* *» » — ^ ,ayersX pnces may be 

entered as a table: 

for 0 > ajq > Zq > bjq corresponding price is Pjq and then 
j=Kq 

Cq(Zq) = 0.5*(Z(|Zq-ajq| + ajq-|Zq-bjq|-bjq) Pjq + 

+ ( | Zq - bKq | - ( Zq - bKq ) )*Pqmax. 
If Zq>0,thenCq(Zq) = 0. 

Kq is number of levels of excavation price at the point Q, 

Pqmax is excavation price for Zq < bKq. (Here Zq, ajq, bjq, bKq < 0). 

If Q is BGP, then 
j=M 

Vq(Zq) =0.5*( H | Zq - ej | + ej - | Zq - fj | - fj )* Fj +■ 
j =1 

+( | Zq - fM | - (Zq - fM))* Fmax), 

7or'o> ej > Zq > fj corresponding price of retaining wall for excavation is Fj. 
M is number of levels of wall for excavation prices. 
Fmax is excavation wall price for Zq < fM <0. 
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If Zq > 0, then Vq(Zq) = 0. 

Wq(Zq) =0.5*(S(|Zq-cj|-cj-|Zq-dj| + dj)* Wj-K|Zq-dL| + (Zq-dL))* 
*Vmax), 

where . 

for 0 < cj < Zq < dj corresponding price of retaining wall for excavation is Wj: 
L is number of levels of fill wall prices. Wmax is wall price for 
Zq > dL. 

If Zq < 0, then Wq(Zq) = 0. 
Finally 

j=N q=N q=N q=N 

D = 0.5*G*G*((I * ( Z Zq + | L Zq | )+0*( 1 1 Zq | - I Zq )), 

q=l q=l q=l Q =1 

I - price of bringing of lack fill from other site to roads system area; 
O - price of removing of redundant fill outside of the roads system area; 
N is total number of GPs on roads system. 
Constraints include five sets: 

a) Constraints of road surface 

Zl = \l*Z'\ + u.l*Z"l -Zle+Xl*Z'e + .ul*Z"e, 
Z2 = \2*Z'2 + u2*Z"2 - Z2e + X2*Z'e + u2*Z"e , 



ZN = \N*Z"N + uN*Z"N - ZNe + XN*Z'e + uN*Z"e , 
Where 

Zq is still decision variable that means relative elevation at Q e GPs for q-l,...N; 

Z'q and Z"q are differences of Z coordinates (proposal - existing) of adjacent CPs 

so that point Q is located between them; 

Zqe is existing elevation at point Q. 

Z'e, Z"e are existing elevations at adjacent CPs. 

Xq, [iq are coefficients of linear interpolation of proposal elevation at point Q. 
Normally there are more than two GPs on each interval between CPs. Therefore 
Z'q and Z"q can be linearly expressed via each pair Zql and Zq2 (except case: 
X. q l=^q2; in this case need use other pair). This allows eliminate Z'q and Z"q from 
surface constraints. Hence elevations of CPs are not used in these constraints 
directly. Instead of them are used their linear expressions via elevations of GPs 
that belong to the interval between two adjacent CPs. However these GPs are all 
located from one side from CP. This obliges introduce additional constraints: 

b) CPs relation constraints: 

Above mentioned linear expressions for Z'q and Z"q via pairs of Zq must be equal 
to such linear expressions for both adjacent intervals of CPs respectively. 

c) Constraints of allowable grades (separately for CPs near to intersection (e.g. 50m) 
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and for other CPs): 

|H' + Z' -(H + Z)|/D< S, 

for all pairs of adjacent CPs; 

where 

H' H existing elevations of pair of adjacent CPs , 

Z\ Z are relative elevations of these CPs; instead Z' and Z implicit their linear 

expressions via corresponding Zq of GPs; 

D is distance along the road between these CPs; 

S is maximal allowable grade for road, or for intersection (e.g. S = 10% or 6%). 

1) Constraints of intersections: (elevations of intersection points of different roads 
must be same). 
H + Z - (H'+ Z') = 0, 

For all CPs at intersections; instead Z' and Z implicit their linear expressions via 
corresponding Zq of GPs; 

Where H, FT are existing elevations at the intersection CPs. 

e) Constraints for lack or excess total fill or excavation volumes: 
q=N q=N 
Z Zq > 0 or Z Zq <0, 
q=l q=l 

The objective function (1.1) contains absolute values. Those absolute values can be 
transformed into linear form. 
In case of 

q=N 

(12) Z Zq > 0, i.e. fill > cut, 
q=l 

the objective (1.1) is reduced as follows: 

q=N j =N 
Z ( Fq(Zq)*G + Wq(Zq) + Cq(Zq)*G + Vq(Zq)) ^-2*G*I * Z Zj. 

q-1 . j=1 

Equation (1 .2) obviously must be joined to constraints. 

Then substitute for each Zq - ai, Zq - bi, Zq - cj, Zq - dj , Zq - ej, Zq - fj 
new variables: 

(*) Zq - ai = Uqail - Uqai2 , i=l,...,Kq, 
(*) Zq-bi = Uqbil -Uqbi2, i=l,...,Kq, 
(*) Zq - bKq = UqbKl - UqbK2; 
(*)Zq-cj = Uqcjl -Uqcj2, j=l,...,L, 
(*) Zq - dj = Uqdjl - Uqdj2, j=l,. . .,L, 
(*) Zq - dL = UqdLl - UqdL2; 
(*) Zq - ej = Uqej 1 - Uqej2, j= 1 , . ,M, 
(*) Zq - fj = Uqfj 1 - Uqfj2, j=l,. . .,M, 
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(*) Zq - flVl = UqfiVll - UqfiVI2. 

In objective instead | Zq -ai | substitute Uqail + Uqai2 etc. So is got LP task because all 
constraints are linear. All constraints also do not include Zq, but only their expressions 
via Uq... l,Uq.--2. 



For example, nonlinear form of Cq(Zq): 
j=Kq 

Cq(Zq) = Z( | Zq - aj | +aj - | Zq - bj | - bj )*Pj + ( | Zq - bKq | - ( Zq - bKq ) )* Pqmax, 
j=l 

will be replaced by linear form: 
j=Kq 

Cq(Zq) = Z(Uqajl+Uqaj2-Uqbjl-Uqbj2 )*Pj+(2*UqbK2)* Pqmax, 
where 

all aj, bj, Pj, bKq, Pqmax are constants and their sums and products do not influence on 
optimization. 

In (*) equations after optimization will be U. . . 1=0 and U. . . 2>=0, 

or U. . .2=0 and U. . . 1>=0 due to objective and implicit constraints make all variables 

>=0. 

The result can be improved by adding CPs and by allowing higher grade. 
Accuracy of the result can be improved by adding density of grid and DTM, 
But this increases dramatically number of variables. 

2. Centerline optimization (Linear programming approach): 

The optimization is reached through linear programming when objective is earthworks 
total cost; variables are differences between proposal roads elevations and existing 
topography elevations at all points of grid on right of way; constraints are present roads 
standards. 

This approach takes on account only CPs i.e. points on centerlines of roads system. 
Grade of road can be changed only at CPs. Recommended to locate CPs at points of 
centerlines where existing profile changes grade. This will improve optimization process. 
On Fig. 3 is presented vertical existing and proposal profiles of road centerline. 
Let Sabo is area of AABO, Sodc is area of AODC. Then 
Sabo- Sodc- AA' x (BA + CD) / 2; 

Where 

AA* is 2d distance between A and D (later Dij); 

AB, DC are signed relative proposal elevations at CPs A and D (later Zij); 
Sabo is fill quantity; 
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Sodc is cut quantity; 

Objective is: difference between cut and fill volumes calculated for roads centerhnes 
only, i.e. for venical profiles: 
j=K i=Nj 

(2 1) Z Z ( Di-lj + Dij)* Zij -> minimum (or max : see below); 
j=l i=l 

where 

K is number of roads; Nj is number of CPs on j-th road of the system; 

Dij is 2d distance between adjacent CPs with relative elevations Zij and Zi+lj on j-th 

road; 

DOj = DNjj= 0 per definition; 

j = 1.....K, 

i=l,...,Nj, 

To satisfy implicit constraints for all variables: >=0 instead relative elevations we used 
Zij+ 10,000. That does not affect on constraints that include only differences of variables, 
but in objective (2. 1) should be added the right hand side member: 
j=K i=Nj 

-Z Z ( Di-lj + Dij)* 10000. 

It should be considered for both cases. When objective is negative it must be maximized. 
When objective is positive it must be minimized. Constraint should be added accordingly: 
j=K i=Nj j=K i=Nj 

Z Z ( Di-lj + Dij)* (Zij +10000) - Z Z ( Di-lj + Dij)* 10000O or >0 
j=l i=l j=l i=l 

Constraints: 

a) Elevations of intersection points of different roads must be same: 
Zljl + Hljl - Zlj2 - Hlj2 = 0, 

Where 

1st point of j 1-th road coincides with 1st point of j2-th road; 
Zljl and Zlj2 are proposal elevations at these points; 
Hljl and Hlj2 are existing elevations at these points. 

b) Grades at all CPs cannot exceed allowable values: 
| Zij + Hij - Zi+lj -Hi+lj | / Dij < Gmax, 

forj = l,...,K;i=l,...,Nj. 



Elevation at each CP should belong to maximal possible allowable interval: 

i=K j=Nj 

ZiOjO < 10000 + max {0, 0.5*(max max ( Hij - HiOjO - DiOjOij * Gmax)}, 

i=l j=l 
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and 

i=K j=Nj 

ZiOjO > 10000 - min {0, 0.5*(min min ( Hij - HiOjO + DiOjOij * Gmax)}, 

i=l j=l 

forjO = 1.....K; i0=l,...,Nj 

Where 

ZiOjO-10000 is proposal relative elevation at the iO-th point on jO-th road; 
DiOjOij is shortest distance between iO-th point on jO-th road to i-th point on j-th road. 
This distance is weighted graph shortest distance. Here weights of graph arcs are 
distances according to road paths. For this purpose we used the Shortest Paths Algorithm 
that described hereafter in paragraph 8. 

3 Centerline optimization (Nonlinear programming approach): 

The optimization is reached through nonlinear programming when objective is 
earthworks total cost; variables are differences between proposal roads elevations and 
existing topography elevations at all points of grid on right of way; constraints are present 
roads standards. 

This approach also takes on account only CPs i.e. points on centerlines of roads system. It 
is recommended locate CPs at points of centerlines where existing profile changes grade. 
This will improve optimization process. , . . , 

Objective is: total sum of unsigned differences between proposal and existing elevations 
calculated for roads centerlines only, i.e. for vertical profiles (Fig. 4): 

A, B, C, D, E are CPs; . . 

AA' is difference between proposal elevation and existing elevation at point A (later Li]). 

j=K i=Nj 

Z Z I Zij | -> minimum; 
j=l i=l 

where 

K is number of roads; Nj is number of CPs on j-th road of the system; 

Dij is distance between adjacent CPs with relative elevations Zij and Zi+lj on j-th road; 

DOj = DNjj= 0 per definition; 

j = l,...,K, 

i = l,...,Nj. 

Constraints: 

a) Elevations of intersection points of different roads must be same: 
Zljl +Hljl-Zlj2-Hlj2 = 0, 
Where: 

1st point of j 1-th road coincides with 1st point of j2-th road; 
Zljl and Zlj2 are proposal elevations at these points; 
Hljl and Hlj2 are existing elevations at these points. 
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b) Maximal grades at all CPs cannot exceed allowable values: 
J Zij + Hif- Zi+lj -Hi+lj | / Dij < Gmax, 
forj = 1,...,K; i=l Nj. 

To transform objective to linear form we use 

(3.1.) Zij = Uijl -Uij2 , 

Uijl,Uij2>=0, 

forj = l,...,K;i=l,...,Kj. 

So in objective instead | Zij | must be Uij 1 + Uij2. 

In (3 . 1 .) equations after optimization will be U. . . 1=0 and U. . . 2>=0, 

or u 2=0 and U. .. 1>=0 due to objective and implicit constraints make all variables 

>=0. 

In all constraints also use for Zij its expression (3.1.). 
4. Refine: 

If CPs are located close each to other (e.g. 10-50m) is recommended to add additional 
constraints: difference between adjacent road grades should not exceed predefined value 
per certain distance (e.g. 0.5% per 10m). Then approximate received optimal centerlines 
by arcs of circles or parabolas or straight lines to satisfy roadway design rules. 

If CPs are located far each to other (e.g. 50-300m) inscribe arcs of circles or parabolas 
between straight intervals according to roadway design rules. 

After the solution obtained, start iteration process of moving entire model along axis Z to 
minimize total cost with taking on account roads cross-sections. The result of 
optimization is set of CPs elevations that define vertical alignment of entire roads system. 

The result can be improved by adding CPs and by allowing higher grade. 
Accuracy of the result can be improved by adding density of grid and DTM. 

5, Fixed points: 

Fixed points means centerline points that must be at definite height. For such points 
would be added fixed points constraints: 
(5.1) H + Z = E, 

where 

Z is relative elevation of the fixed point; 
H is existing elevation at this point; 
E is required height at the fixed point. 

Equation (5.1) take place for each fixed point. It means that each fixed point must be 
defined as CP. If there are two fixed points or more, in each optimization method, it is 
necessary to ensure, before optimization, that grades between the fixed points do not 
exceed the maximal allowable value. This checking must be done for the shortest 
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distance between fixed points. For that purpose we used the Shortest Paths Algorithm that 
described hereafter in paragraph 8. 

6. Subdivision: 

If the optimization process includes earthworks of subdivision lots then it is necessary 
added to the model the proposal sections and boundaries of lots. 

7. Optimization Algorithm for excavation material Transportation: 

After optimization of vertical profiles of roads is done can be received tables of fill-and- 
cut quantities according to predefined intervals along roads (say, each 50m). So the work 
program for transportation of excavation material at preliminary planning stage should be 
created. The invention allows receiving optimal work program. In our model the 
suppliers are sites with excess of excavation material and demanders are sites with lack of 
excavation material. The objective is minimum of sum of 

excavation weights multiplied by distances: 
i=Nj=M 

Z Z Dij * Xij -> minimum; 
i=l j=l 

where 

Xij is the decision variable that means fill quantity will be moved from i-th supplier to j- 
th demander; i= 1 , . . . ,N, j = 1 , . . . ,M; 

where N is number of suppliers and M is number of demanders; 

Dij is a shortest distance from i-th supplier to j-th demander; Dij are input data for this 
optimization process. 

Constraints are: 

a) Move not more than each demander needs: 
i=N 

ZXij<Pj;j=l,...,M; 
i=l 

Pj is fill of j-th demander; 

b) Move not more than each supplier has: 
j=M 

ZXij<Qi; i=l,...,N; 

Qi is fill of i-th supplier; 

c) But move all that possible: 
i=N j=M j=M j=N 

Z Z Xij = min ( Z Pj , Z Qi ). 
i=l j=l j=l i=l 
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The shortest distance between two points is defined through shortest paths 
algorithm described in paragraph 8. 

8 The Shortest Paths Algorithm: 

In this invention we used the shortest paths algorithm for Centerhne optimization -Linear 
programming approach (Paragraph 2), fixed points (Paragraph 5) and transportation 
(Paragraph 7). 

DESCRIPTION OF THE DRAWING: 

Fig. 1 : Scheme of the optimization process. 

Relation between grid and control points for grid optimization. 
Fill - Cut geometrical presentation for centerline linear model objective function. 
Relative proposal unsigned elevations for centerline nonlinear model objective 
function. 



Fio 2 
Fig. 3 
Fis. 4 



CLAIMS: 

1 Using process of optimization to obtain the optimal: points of intersection, 

vertical alignment, cross-sections, quantity, cost estimation of the entire project 
and display result, according to subject plan, topography, roads standards and 
units cost (include different prices for different geologic layers) in roads 
planning. 

2. The process of claim 1 wherein "optimization" comprises linear, nonlinear 
programming. 

3 The process of claim 1 wherein Optimal means an objective function that 
minimizes total cost depends on a set of variables while satisfying some 
constraints. 

4. The process of claim 1 wherein display result comprises 2D & 3D CAD or 
alphanumeric information. 

5. The process of claim 1 wherein Topography comprises DTM, TIN 

6. The process of claim 1 wherein Road comprises highway, railway, pedestrian. 

7. The process of claim 1 wherein comprises subdivision. 

8. Using geologic layers data for the optimization mentioned in claim 1 . 
9 Using process of optimization for the entire project. 

10. Using process shortest paths algorithm for roads system optimization. 



13 

1 1. The process of claim 10 wherein comprises roads system with fix points. 

12. Using process of optimization to obtain optimal excavation material 
transportation. 

13. Utilizes the information generated by a computer-aided design (CAD) program as 

input to the optimization algorithm for road design. ♦ 

14. Using process of comparison different optimal land use layouts to select the one / 
that minimizes total cost for road design purpose. 
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